A HYBRID OF DARBOUX'S METHOD 
AND SINGULARITY ANALYSIS IN COMBINATORIAL 

ASYMPTOTICS 



PHILIPPE FLAJOLET, ERIC FUSY, XAVIER GOURDON, DANIEL PANARIO, AND 

NICOLAS POUYANNE 



Abstract. A "hybrid method" , dedicated to asymptotic coefficient extraction 
in combinatorial generating functions, is presented, which combines Darboux's 
method and singularity analysis theory. This hybrid method applies to func- 
tions that remain of moderate growth near the unit circle and satisfy suitable 
smoothness assumptions — this, even in the case when the unit circle is a natu- 
ral boundary. A prime application is to coefficients of several types of infinite 
product generating functions, for which full asymptotic expansions (involving 
periodic fluctuations at higher orders) can be derived. Examples relative to 
permutations, trees, and polynomials over finite fields are treated in this way. 



Introduction 

A few enumerative problems of combinatorial theory lead to generating functions 
that are expressed as infinite products and admit the unit circle as a natural bound- 
ary. Functions with a fast growth near the unit circle are usually amenable to the 
saddle point method, a famous example being the integer partition generating func- 
tion. We consider here functions of moderate growth, which are outside the scope 
of the saddle point method. We do so in the case where neither singularity analysis 
nor Darboux's method is directly applicable, but the function to be analysed can 
be factored into the product of an elementary function with isolated singularities 
and a sufficiently smooth factor on the unit circle. Such decompositions are often 
attached to infinite products exhibiting a regular enough structure and are easily 
obtained by the introduction of suitable convergence factors. Under such condi- 
tions, we prove that coefficients admit full asymptotic expansions involving powers 
of logarithms and descending powers of the index n, as well as periodically vary- 
ing coefficients. Applications are given to the following combinatorial-probabilistic 
problems: the enumeration of permutations with distinct cycle lengths, the prob- 
ability that two permutations have the same cycle-length profile, the number of 
permutations admitting an mth root, the probability that a polynomial over a fi- 
nite field has factors of distinct degrees, and the number of forests composed of 
trees of different sizes. 

Plan of the paper. We start by recalling in Sectionnthe principles of two clas- 
sical methods dedicated to coefficient extraction in combinatorial generating func- 
tions, namely Darboux's method and singularity analysis, which are central to our 
subsequent developments. The hybrid method per se forms the subject of Section|21 
where our main result. Theorem El is established. Section 01 treats the asymptotic 
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enumeration of permutations having distinct cycle sizes: this serves to illustrate 
in detail the hybrid method at work. Section 01 discusses more succinctly further 
combinatorial problems leading to generating functions with a natural boundary — 
these are relative to permutations, forests, and polynomials over finite fields. A 
brief perspective is offered in our concluding section, Section [S] 

1. DARBOUX'S method and singularity ANALYSIS 

In this section, we gather some previously known facts about Darboux's method, 
singularity analysis, and basic properties of analytic functions that are central to 
our subsequent analyses. 

1.1. Functions of finite order. Throughout this study, we consider analytic func- 
tions whose expansion at the origin has a finite radius of convergence, that is, func- 
tions with singularities at a finite distance from the origin. By a simple scaling of 
the independent variable, we may restrict attention to function that are analytic 
in the open unit disc D but not in the closed unit disc D. What our analysis a 
priori excludes are thus: (i) entire functions; (ii) purely divergent series. (For such 
excluded cases, either the saddle point method or ad hoc manipulations of divergent 
series are often instrumental in gaining access to coefficients |2|[Tll|2n|-) Further- 
more we restrict attention to functions that remain of moderate growth near the 
unit circle in the following sense. 

Definition 1. A function f{z) analytic in the open unit disc D is said to be of 
global order a < if 

f{z)^o{{i-\z\r) (|z|<i), 

that is, there exists an absolute constant M such that \f{z)\ < M(l — |z|)° for all z 
satisfying \z\ < 1. 

This definition typically excludes the partition generating function 

oo 

- n 

k=l 

which is of infinite order and to which the saddle point method (as well as a good 
deal more) is applicable ^I^IEOI- In contrast, a function like 



is of global order a — —i, while 



\ - z'^ 



are of global order a — 0. 

We observe, though we do not make use of the fact, that a function f{z) of global 
order a < has coefficients that satisfy [z'^]f{z) = 0(n~°). The proof results from 
trivial bounds applied to Cauchy's integral form 

upon integrating along the contour C: \z\ — 1 — n^^. (In \7\, Braaksma and Stark 
present an interesting discussion leading to refined estimates of the 0{n~°-) bound.) 
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1.2. Log-power functions. What we address here is the asymptotic analysis of 
functions whose local behaviour at designated points involves a combination of 
logarithms and powers (of possibly fractional exponent). For the sake of notational 
simplicity, we write 

L{z) log Y^T^- 
Simplifying the theory to what is needed here, we set: 
Definition 2. A log-power function at 1 is a finite sum of the form 

r 

o{z) = Y^c,{Uz)){i-zr\ 

k=l 

where ai < ■ ■ ■ < and each cj, is a polynomial. A log-power function at a finite 
set of points Z = {Ci, . . . , Cm}; a finite sum 

m / 

where each aj is a log-power function at 1. 

In what follows, we shall only need to consider the case where the Q lie on the unit 
disc: ICjl — 1. 

It has been known for a long time (see, e.g., Jungen's 1931 paper, ref. |22| . 
and ^1 for wide extensions) that the coefficient of index n in a log-power 
function admits a full asymptotic expansion in descending powers of n. 

Lemma 1 (Coefficients of log-powers). The expansion of the coefficient of a log- 
power function is computable by the two rules: 

1 (—a) I y—a) 

(2) [z-][l-zrUzf - (-l)'=^([z"](l-z)") 



r^-a-i Q(a + l)n-"-2 



_r(-a) r(-a) 

The general shape of the expansion is thus 

[z"](l-z)"L(^)'= ^ .^^n-"-i(logn)'= (a ^ Z>o) 

[z"](l-z)''L(z)'= - (-l)''fc(r!)n-'^-i(logn)'=-i (r G Z>o, fc G Z>i). 

n — *+oo 

In the last case, the term involving (logn)*"' disappears as its coefficient is l/r(— r) = 
0. In essence, smaller functions at a singularity have asymptotically smaller coef- 
ficients and logarithmic factors in a function are reflected by logarithmic terms in 
the coefficients' expansion; for instance, 

_L{z)_ logn-t-7 + 21og2 logn + 7 + 21og2 
[z"](l-z)L(z)2 ~ ^(log7i + 7-l)-— (21og7i + 27-5)-h--- . 



When supplemented by the rule 
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Lemma ^ makes it effectively possible to determine the asymptotic behaviour of 
coefficients of all log-power functions. 

1.3. Smooth functions and Darboux's method. Once the coefficients of func- 
tions in some basic scale are known, there remains to translate error terms. Pre- 
cisely, we consider in this article functions of the form 

/(^) = I](z)-hi?(z), 

and need conditions that enable us to estimate the coefficients of the error term 
R{z). Two conditions are classically available: one based on smoothness (i.e., 
differentiability) is summarized here, following classical authors (e.g., the 
other based on growth conditions and analytic continuation is discussed in the 
next subsection. 

Definition 3. Let h{z) be analytic in \z\ < 1 and s be a nonnegative integer. The 
function h{z) is said to be C-smooth^ on the unit disc ( or of class ) if, for 
all k = . . s, its kth derivative h''^\z) defined for \z\ < 1 admits a continuous 
extension on \z\ < 1. 

For instance, a function of the form 

h{z) = /i„z" with /i„ = Oin-'-^-^), 

n>0 

for some (5 > and s e Z>o, is C-smooth (both in the standard sense and in the 
strong sense). Conversely, the fact that smoother functions have asymptotically 
smaller coefficients lies at the heart of Darboux's method. 

Lemma 2 (Darboux's transfer). If h{z) is -smooth, then 

[z"]/i(z) = o{n-'). 

Proof. One has, by Cauchy's coefRcient formula and continuity of h{z): 

/i(e*^)e-'"^ de. 

When s = 0, the statement results directly from the Riemann-Lebesgue theo- 
rem [331 P- 109]. When s > 0, the estimate results from s successive integra- 
tions by parts followed by the Riemann-Lebesgue argument. See Olver's book 
p. 309-310] for a neat discussion. □ 

Definition 4. A function Q{z) analytic in the open unit disc D is said to admit a 
log-power expansion of class C* if there exist a finite set of points Z = {^i, . . . , t^m} 
on the unit circle \z\ — 1 and a log-power function E(z) at the set of points Z such 
that Q{z) — S(z) is C* -smooth on the unit circle. 

Lemma 3 (Darboux's method). If Q{z) admits a log-power expansion of class C* 
with S(z) an associated log-power function, its coefficients satisfy 

[z"]Q(z) = [z"]E(z)-fo(n-*). 

function h{z) is said to be weakly C" -smooth if it admits a continuous extension to the 
closed unit disc \z\ < 1 and tiie function g{9) := h(e**) is s times continuously differentiable. This 
weaker notion suffices for Lemmas and Theorem [Tl 
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Proof. One has Q — II + R, with R being C* smooth. The coefficients of R are 
estimated by Lemma El □ 

Consider for instance 

y'l — z vl — z 

Both are of global order — 5 in the sense of Definition ^ By making use of the 
analytic expansion of at 1, one finds 



[-^^-eVT^)+R,i.), 



where Ri{z), which is of the order of (1 — z)^^^ as z ^ 1^, is C^-smooth. The sum 
of the first two terms (in parentheses) constitutes S(z), in this case with Z = {1}. 
Similarly, for Q2{z), by making use of expansions at the elements of Z = { — 1, +1}, 
one finds 

Q2{z) = (^fi= - ^V2VT~z + -^VT+i] + i?2(z), 
yy/l~z 4 eV2 J 

where R2{z) is also C^'^-smooth. Accordingly, we find: 



n I \ T^n \ n 



The next term in the asymptotic expansion of [z"]Q2 involves a linear combination 
of n~'^/^ and (— where the latter term reflects the singularity at z = — 1. 
Such calculations are typical of what we shall encounter later. 

1.4. Singularity analysis. What we refer to as singularity analysis is a technology 
developed by Flajolet and Odlyzko |14l \'M)\. with further additions to be found 
in jlOl 111! I15j . It applies to a function with a finite number of singularities on the 
boundary of its disc of convergence. Our description closely follows Chapter VI of 
the latest edition of Analytic Combinatorics |15) . 

Singularity analysis theory adds to Lemma^the theorem that, under conditions 
of analytic continuation, O- and o-error terms can be similarly transferred to coef- 
ficients. Define a A-domain associated to two parameters R > 1 (the radius) and 
0e (0, f) (the angle) by 

A(i?, (f>) := {z I \z\< R, (j)< arg(z ~ 1) < 2tt - (j), z ^ l} 

where arg(w) denotes the argument of w taken here in the interval [0, 27r[. By 
definition a A-domain properly contains the unit disc, since </> < ^. (Details of the 
values of R, cj) are immaterial as long as i? > 1 and cj) < 

The following definition is in a way the counterpart of smoothness (Definition^ 
for singularity analysis of functions with isolated singularities. 

Definition 5. Let h(z) be analytic in \z\ < 1 and have isolated singularities on 
the unit circle at Z — {^1, . . . , Cm}- Let t be a real number. The function h{z) is 
said to admit a log-power expansion of type O* ( relative to Z) if the following two 
conditions are satisfied: 

— The function h{z) is analytically continuable to an indented domain !D = 
Pli'lilCj ■ ^); ^ some A-domain. 
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— There exists a log-power function ^{z) :— X^JLi ^j{^Kj) such that, for each 
Cj e Z, one has 

(4) hiz)-a,iz/Q) = 0{{z-C,Y), 
as z C,j in (Q ■ A). 

Observe that "^iz) is a priori uniquely determined only up to 0{{z — Cj)*) terms. 
The minimal function (with respect to the number of monomials) satisfying Q is 
called the singular part of h(z) (up to O* terms). 

A basic result of singularity analysis theory enables us to extract coefficients of 
functions that admit of such expansions. 

Lemma 4 (Singularity analysis method). Let Z = be a finite set 

of points on the unit circle, and let P{z) he a function that admits a log-power 
expansion of type O* relative to Z , with singular part ^(z). Then, the coefficients 
of h satisfy 

(5) [z"]P(z) = [z"]S(^)+0(n-*-i). 

Proof. The proof of Lemma^jstarts from Cauchy's integral formula 10 and makes 
use of the contour C that lies at distance — of the boundary of the analyticity 
domain, 2) = fljli (0 - A). See [H [El for details. □ 

1.5. Polylogarithms. For future reference (see especially Section OJ, we gather 
here facts relative to the polylogarithm function Lii,(z), which is defined for any 

OO 

(6) Li.W:=E^- 
One has in particular 



n— 1 



Z 1 

Lio(z) = , Lii(z)=log- =L(z). 

1 — z 1 — z 

In the most basic applications, one encounters polylogarithms of integer index, but 
in this paper (see the example of dissimilar forests in Section ^ , the more general 
case of a real index v is also needed. 

Lemma 5 (Singularities of polylogarithms). For any index v d C, the polyloga- 
rithm Li^{z) is analytically continuable to the slit plane C \ K>i. If v ^ in E 
the singular expansion of Li™ (z) near the singularity z — I is given by 



(7) 



Li™W = / \. T--Hlogr-g.„-i)+ V ^—r^Cim-jy 

TO — 1)! '^-^ V. 



e=i 

For V not an integer, the singular expansion of\A^{z) is 

(8) Li.(z) ~ r(i - i.)T^-' + ^ ^^C(^ - J>^' 
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The representations are given as a composition of two explicit series. The expan- 
sions involve both the harmonic number Hm and the Riemann zeta function (^(s) 
defined by 

11 1 ,,111 

(C(s), originally defined in the half-plan 3ff(s) > 1, is analytically continuable to 
C \ {1} by virtue of its classical functional equation). 



Proof. First in the case of an integer index m g Z>2, since Lim(^) is an iterated 
integral of Lii(z), it is analytically continuable to the complex plane slit along 
the ray [l,-|-oo[. By this device, its expansion at the singularity z — I can be 
determined, resulting in Q. (The representation in Q is in fact exact and not 
merely asymptotic. It has been obtained by Zagier and Cohen in |27[ p. 387], and 
is known to the symbolic manipulation system Maple.) 

For V not an integer, analytic continuation derives from a Lindelof integral repre- 
sentation discussed by Ford in ^16j . The singular expansion, valid in the slit plane, 
was established in to which we refer for details. □ 

In the sequel, we also make use of smoothness properties of polylogarithms. 
Clearly, Lifc(z) is C'^~^-smooth in the sense of Definition |21 A simple computation 
of coefficients shows that any sum 

Skiz) = J2rii) [Li,(z^)-Li,(l)] 

e>k 

with r{x) polynomially bounded in x, is C'^^^-smooth. Many similar sums are 
encountered later, starting with those in Equations H25(l and H26|l . 

2. The hybrid method 

The heart of the matter is the treatment of functions analytic in the open unit 
disc that can, at least partially, be "de-singularized" by means of log-power func- 
tions. 

2.1. Basic technology. Our first theorem, which essentially relies on the Darboux 
technology, serves as a stepping stone towards the proof of our main statement. 
Theorem El below. 

Theorem 1. Let f{z) be analytic in the open unit disc D, of global order a < 0, and 
such that it admits a factorization f = P ■ Q, with P, Q analytic in D. Assume the 
following conditions on P and Q, relative to a finite set of points Z = {(^i, . . . , Cm} 
on the unit circle: 

Ci: The "Darboux factor" Q{z) is C -smooth on the unit circle (s G Z>Oy). 

C2: The "singular factor" P{z) admits, for some nonnegative integer t, a log- 
power expansion relative to Z , P = P + R (with P the log-power function 
and R the smooth term), that is of class C*. 

Assume also the inequality (with [xj the integer part function): 

C3: t > uq, where 

(9) .„ - ^ ' L«J 
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Let Co ~ ^ ^""^ j . If H denotes the Hermite interpolation polynomial such that all 
its derivatives of order 0, . . . , Cq — 1 coincide with those of Q at each of the points 
Ci, • ■ ■ ,Cm, one has 

(10) = n (P(z) • H{z)) + o(n-«''). 

Since P{z)-H{z) is itself a log-power function, the asymptotic form of its coefficients 
is explicitly provided by Lemma ^ 

Proof. Let c < s be a positive integer whose precise value will be adjusted at the 
end of the proof. First, we decompose Q as 

Q = Q + S, 

where Q is the polynomial of minimal degree such that all its derivatives of order 
0, . . . , c — 1 at each of the points Ci, • ■ ■ , Cm coincide with those of Q: 



(11) 



, < i < c, 1 < J < 

The classical process of Hermite interpolation [21] produces such a polynomial, 
whose degree is at most cm — 1. Since Q is C°°-smooth, the quantity S = Q — Q 
is C*-smooth. This function S is also "flat" , in the sense that it has a contact of 
high order with at each of the points Q . 

We now operate with the decomposition 

(12) f = P-Q + P-S + R-Q + R-S, 
and proceed to examine the coefficient of z" in each term. 

— The product P ■ Q. Since P is a log-power function and Q a polynomial, 
the coefficient of z" in the product admits, by Lemma ^ a- complete descending 
expansion with terms in the scale {ri^''(logr7,)'^}, which we write concisely as 

(13) [z"]P -Q e {n-^^ {log n)^ \ fc G Z>o, /3 G M} . 

— The product P ■ S. This is where the Hermite interpolation polynomial Q 
plays its part. From the construction of Q, there results that S — Q ~ Q has all 
its derivatives of order 0, . . . , c — 1 vanishing at each of the points Ci, • ■ • , Cm- This 
guarantees the existence of a factorization 

m 

S{z) = Q{z) - Q{z) = k{z) l[{z~ 0)^ 

where k{z) is now C^^'^-smooth (division decreases the degree of smoothness). 
Then, in the factorization 

P-S^ h-f[{z-Qr] -niz), 



Hermite interpolation extends the usual process of Lagrange interpolation, by allowing for 
higher contact between a function and its interpolating polynomial at a designated set of points. 
A lucid construction is found in Hildebrand's treatise .21. §8.2], 
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the quantity P is, near a of order at most 0{z — (with a the global order of 
/). Thus, PS/k is at least C"-smooth, with v ■.= [c + a\. Since • C« C ^•^'"(P'?), 
Darboux's method (Lemma (JJ yields 

(14) [z"]P • = o (^n-"^'^)) , u{c) ■.= mm{[c + a\,s-c). 

— The product R ■ Q. This quantity is of class C* by elementary product rules. 
Hence, by Darboux's method, 

(15) -Q = 0(71-*) . 

— The product R S. This product is of class and, by Darboux's method 
again, 

(16) •S' = o(n- . 

It now only remains to collect the effect of the various error terms of (jl4|l , (|15|l , 
and Hlf)|l in the decomposition 1)12(1 : 

- min(s,i) \ 



']/ = ([z"]P • q) + o(n-"(=)) + o{n'') + o{n- 



Given the condition t > uq in C3, the last two terms are o(n "°). A choice, which 
maximizes u{c) (as defined in ((14() ') and suffices for our purposes, is 



(17) 



s- [aj 



corresponding to u{cq) = 



[a\ 



Uq. 



The statement then results from the choice of c = cq, as well as uq = u{co) and 
H(z) := Q{z), the corresponding Hermite interpolation polynomial. □ 

2.2. Hybridization. Theorem Q] is largely to be regarded as an existence result: 
due to the factorization and the presence of a Hermite interpolation polynomial, it 
is not well suited for effectively deriving asymptotic expansions. In this subsection, 
we develop the hybrid method per se, which makes it possible to operate directly 
with a small number of radial expansions of the function whose coefficients are to 
be estimated. 

Definition 6. Let f{z) be analytic in the open unit disc. For ( a point on the unit 
circle, we define the radial expansion of f at Q with order t G M as the smallest ( in 
terms of the number of monomials) log-power function a{z) at C, provided it exists, 
.such that 

f{z)^aiz)+0{iz-CY), 
when z = (1 — x)C and x tends to 0+. The quantity cr(z) is written 

asymp(/(2;),C,i). 

The interest of radial expansions is to a large extent a computational one, as 
these are often accessible via common methods of asymptotic analysis while vari- 
ous series rearrangements from within the unit circle are granted by analyticity. In 
contrast, the task of estimating directly a function f{z) as z ^ C on the unit circle 
may be technically more demanding. Our main theorem is accordingly expressed 
in terms of such radial expansions and, after the necessary conditions on the gener- 
ating function have been verified, it provides an algorithm (Equation H18|l ) for the 
determination of the asymptotic form of coefficients. 
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Theorem 2 (Hybrid method). Let f{z) he analytic in the open unit disc D, of 
finite global order a < 0, and such that it admits a factorization f = P ■ Q , with 
P,Q analytic in D. Assume the following conditions on P and Q, relative to a 
finite set of points Z = {Ci, . . . , (^m} on the unit circle: 

Di: The "Darboux factor" Q{z) is C -smooth on the unit circle (s G Z>o/ 
D2: The "singular factor" P{z) is analytically continuable to an indented do- 
main of the form J) = fXjLi iCj ' ^)- For some non-negative real number 
to, it admits, at any Q G Z , an asymptotic expansion of the form 

P{z) = a.iz/Q,) + O ((z - 0)*") {z -> 0, z e S), 

where (Tj{z) is a log-power Junction at 1. 

Assume also the inequality: 

D3: <o> L^J. 

The 



Then f admits a radial expansion at any Q G Z with order uo 
coefficients of f{z) satisfy: 

[2"]/(2) = [z"]A(z) + o(n-"«), 

(18) " 

where A(z) := ^ asynip(/(z), Cj, wq). 



2 



Proof. Let us denote by 



S(z) = > a,(z/C,) 



J'' 
1 



the sum of the singular parts of P at the points of Q. The difference R :— P — Y.\s 
C*-smooth for any integer t satisfying t < to (in particular, we can choose t — uq, 
this by assumption D2 . The singular factor P has thus been re-expressed as the sum 
of a singular part E and a smooth part R. The conditions of Theorem ^a-re then 
precisely satisfied by the product PQ, the inequality D3 implying condition C3, so 
that one has by (|10() 

(19) n/(z) = [z"]S(z)ff(z)+o(n-"«), 

where H is the Hermite polynomial associated with Q that is described in the proof 
of Theorem and uq is given by lO ■ 

In order to complete the proof, there remains to verify that, in the coefficient 
extraction process of (|19|) above, the quantity EiJ can be replaced by A(z). 

We have 

(20) [zimH{z) - ^[z"]a,(z/0)i/(^). 

j 

Now, near each Q, we have (with cq = |^^— 5^ according to ^17^ ) 

a,{z/Q) = P(z) + O((z-0)*°) 

(21) H{z) = Q(z) + O((z-0r«) 
P{z) = 0{{z~Qr), 
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respectively by assumption D2, by the high order contact of H with Q due to the 
Hermite interpolation construction, and by the global order property of f{z). There 
results from Equation H21|) . condition D3, and the value of cq in (|17|l that 

a,{z/QH{z) = asymp(P(z)g(z),0,?^o) + 0{[z - QD, 
The proof, given 1)19(1 and ((2()|l . is now complete. □ 

Thanks to Theorem |21 in order to analyse the coefficients of a function /, the 
following two steps are sufRcient. 

(i) Establish the existence of a proper factorization / = P Q. Usually, a crude 

analysis is sufRcient for this purpose. 
(m) Analyse separately the asymptotic character of f{z) as z tends radially to 
a few distinguished points, those of Z . 

As asserted by Theorem |3 it then becomes possible to proceed with the analysis 
of the coefficients as though the function / satisfied the conditions of 

singularity analysis (whereas in general f{z) admits the unit circle as a natural 
boundary). 

Manstavicius [28| develops an alternative approach that requires conditions on 
generating functions in the disc of convergence, but only weak smoothness on the 
circumference. His results are however not clearly adapted to deriving symptotic 
expansions beyond the main terms. 

3. Permutations with distinct cycle sizes 

The function 

fe=i ^ ^ 

has been studied by Greene and Knuth 19 , in relation to a problem relative to 
factorization of polynomials over finite fields that we treat later. As is readily rec- 
ognized from first principles of combinatorial analysis 15, 17, 36, 40 , the coefficient 
[z"]/(z) represents the probability that, in a random permutation of size n, all cycle 
lengths are distinct. One has 

z'^ Z^ 

/W = l + - + ^ + 5^ + 14- + 74- + ..., 

and the coefficients constitute the sequence EIS:A007838^ In [El §4.1.6], the 

authors devote some seven pages (pp. 52-58) to the derivation of the estimate (7 
is Euler's constant) 

(22) f^:=[,^]f(^) = e-^ + — + o(^-^], 

starting with a Tauberian argument and repeatedly using bootstrapping. In our 
treatment below, we recycle some of their calculations, though our asymptotic 
technology is fundamentally different. 



We shall use the notation EIS:xxxxxx to represent a sequence indexed in the Encyclopedia 
of Integer Sequences |34| . 
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Global order. The first task in our perspective is to determine the global order 
of f{z). The following chain of calculations, 



/(.) = n^^^^'^nfi+y) 

I — 1 I — 1 V / 




-z^/k 










k) 


k>l 


^3k 
fc3 



fc=l 

(23) = T^^^P 

1 

= cxp 

shows f{z) to be of global order —1. It is based on the usual introduction of 
convergence factors, the exp-log transformation {X = exp(logX)), and finally the 
logarithmic expansion. 

Note that this preliminary determination of global order only gives the useless 
bound /„ = 0{n). Actually, from the infinite product expression of the Gamma 
function [21] (or from a direct calculation, as in JHI)) there results that 



(24) = n fi + ^ 



k>l 



hence, from the second line of (|23|l . 

/(^) - , 

z^l- 1 — Z 

which is compatible with (|22|l . but far from sufficient to imply it. 

The hybrid method. Given the last line of H23() . which is re-expressed in terms 
of polylogarithms as 

(25) /(z) = cxp f-i Li2(z2) + 1 Li3(z3) _ . 

1 — z V 2 3 

the right factorization of /(z) is obtained transparently. Define 

(26) U{z):= J2 ^^7^Li^(^'), Viz) := tl^llueiz'), 

l<f<s+l s+2<l 

SO that 

(27) /(z) = e^(") -e^f"). 

Clearly V{z) is C^-smooth, and so is e^*^^^ given usual rules of differentiation. Thus 
Q := is our Darboux factor. The first factor P := eP satisfies the condition 
of Theorem |21 it is the singular factor and it can be expanded to any order t 
of smallncss. Consequently, the hybrid method is applicable and can provide an 
asymptotic expansion of [z"]/(z) to any predetermined degree of accuracy. 

The nature of the full expansion. Given the existence of factorizations of type 1)27(1 
with an arbitrary degree of smoothness (for V) and smallness (for U), it is pos- 
sible to organize the calculations as follows: take the primitive roots of unity in 
sequence, for orders 1,2,3,.... Given such a root rj of order £, each radial restric- 
tion admits a full asymptotic expansion in descending powers of (1 — z/rj) tempered 
by polynomials in log(l — z/77). Such an expansion can be translated formally into 
a full expansion in powers of tempered by polynomials in logn and multiplied 
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by 77 All the terms collected in this way are bound to occur in the asymptotic 
expansion of /„ = [z"]/(z). 

For the sequel, it proves convenient first to adjust the expansion \2'A\i . by taking 
out the (1 + z) factor. We find, by the same techniques as in H23|l 



(28) f{z) - e-^-— exp V [Li^Cz") - 

1 — z \ ^ — ' m 

\m>2 



and start with the expansion as z — > 1, then consider in turn z = —1, z — u^uP' 
(with Lo — exp(2i7r/3)), and finally z = ry, a primitive £th root of unity. 

The expansion at z = \ . Calculations simplify a bit if we set 

z — e^'', T = — log z, 

as in |(7J). By summing the singular expansions of polylogarithms Q, one arrives 
at an asymptotic expansion as r ^ 0+ of the form: 

(29) /(e--) = e-^'^ exp (-a(T) log r + /3(t) + <5i(t) + 52{t) - e(T)) . 

1 — e ^ 

There, the first two terms in the exponential correspond to summing the special 
terms in the singular expansions of polylogarithms Q: 

<^i^)-=Y.- n ' /3M=zI^ n (i/„~iog(m + i)). 

^ — ' m! ^ — ' to! 

m>l m>l 



The last three terms inside the exponential of (|29|l arise from summation over values 
of TO > 2 of the regular part of , namely, 

E^F^ E ^c(™-.)(-)^ 

m>2 j>0, j^m-l •' 

upon distinguishing between the three cases: j < to — 1 (giving rise to Si{t)), 
j = m — 1 (giving e(T)), and j > to — 1 (giving S2{t)), and exchanging the order of 
summations. The calculation of e(T) and (5i(t) is immediate. First 

(30) 5i(r) 5,,, = ^(-l)™-W-i(C(TO-j)-l), 

j>2 m=2 

and since ^ values at negative integers are rational numbers, the expansion of <5i(t) 
involves only rational coefficients. Next, one finds 



e(r)=5](TO + l)"-i — 
to! 



a variant of the Lambert and Cayley functions. Finally, the function (52 (t) has 
coefficients a priori given by sums like in (|3()|l , but with the summation extending 
to TO > + 2: 

3>0 m>]+2 



14 



FLAJOLET, FUSY, GOURDON, PANARIO, AND POUYANNE 



Each infinite sum in the expansion of 82 is expressibie in finite form: it suffices 
to start from the known expansion of + s) at s = 0, which gives ("(/"(s) is the 
logarithmic derivative of the Gamma function) 

(31) ^(1 + ,) + 7 - = (C(2) - l)s - (C(3) - l)s^ + • • • , 

1 + s 

and differentiate an arbitrary number of times with respect to s, then finally set 
s = 1. One finds for instance, in this way, 

-52,0 --7 -log 2 + 1, '52,2 - 2C(3) - y + I 

Equation 129|) provides a complete algorithm for expanding /(z) as z ^ 1^. The 
first few terms found are 



f{z) = e-V- Iog(l-z)-Iog2 + -(l-z)log"(l-z) 

\ 1 — z 2 

+ (Iog 2 - 2)(1 - z) log(l - z) + 0(1 - z)bigg). 

Then, an application of Theorem |21 yields the terms in the asymptotic expansion of 
/„ arising from the singularity z = 1: 

(32) jf^ =e-' + — + ^{-\ogn + cz.o) + ^ (log^ n + 03,1 log n + 03,0) + ■■■ 
where 

.33. C2,o = -1 -7 + log2, C3,i =4 + 27-21og2, 

C3,o = 1 + 47 - log2 - 31og3 + log^ 2 - ^ + 72 - 27log2. 

From preceding considerations, the coefficients all lie in the ring generated by 
log2,log3, . . ., 7, TT, and C(3),C(5), . . .. 

The terms given in (|32|l provide quite a good approximation. Figure ^ displays 

^11 

a comparison between /„ and its asymptotic approximation /A , up to terms of 
order n~^. We find 

|i?„| <22, 

for all 1 < ?i < 1000, and expect the bound to remain valid for all rt > 1. 



= + \Rn\ 



Expansions at z — —\ and at z = uj^lo"^ . We shall content ourselves with brief in- 
dications on the shape of the corresponding singular expansions. Note that Figure^ 
clearly indicates the presence of a term of the form 

VL{n) 

in the full asymptotic expansion of /„, where J7(n) is a function with period 6. This 
motivates an investigation of the behaviour of /(z) near primitive square and cube 
roots of unity. 

Start with z —>■ — 1+. The definition of f[z) implies that 

/(.)_-^^(i+.)n(i+^^) 

k=2 ^ ^ 
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Figure 1. Permutations with distinct cycle lengths: the approx- 
imation error as measured by n^ifn/fn^ — 1), with fll^ truncated 
after n'^ terms, for n = 1 . . 200 (left) and for n = 180 . . 200 (right). 



(the infinite product telescopes). Set Z = 1 + z and restrict attention to the type 
of the expansions at z = — 1. Only half of the polylogarithms in H28|) are singular, 
so that the expansion at z = —1 is of the type 

Zexp{ZlogZ + Z + + Z^logZ + ■ ■ ■) ^ Z + Z^ logZ + Z^ + Z'^ log^ Z + ■ ■ ■ . 

There, we have replaced all unspecified coefficients by the constant 1 for readability. 
This singular form results in a contribution to the asymptotic form of /„: 



J[-l] ^ ^_^yn f^j^ d4,llogn + rf4,0 



(Compared to roots of unity of higher order, the case z = — 1 is special, because of 
the factor (1 -I- z) explicitly present in the definition of /(z).) A simple calculation 
shows that da = 2, so that 

Next, let ui = e^*^/'^ and set Z := (1 — z/uj). The type of the expansion at z = w 

is 

/(z) ^ /(c^)exp(Z + Z2iogZ + z3 + ...) 

- f{uj) {1 + Z + Z'^logZ + Z^ + ZHogZ + Z'^log^ Z + ■ ■ ■) , 

since now every third poly logarithm is singular at z = lu. This induces a contribu- 
tion of the form 

fluj] _ , -n f 64,1 l0gn + 64,0 , 
In ^ ^ ^ H 4 r ■ 

arising from z = uj, and similarly, for a conjugate contribution arising from w^. 
Another simple calculation shows that 

63 - 3/(l^), 

and leaves us with the task of estimating /(w). The use of the formula, 

r(i + a)6-^'^ 



n 1 



k>l 



r(l + s + a) ' 
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Figure 2. Permutations with distinct cycle lengths: the approx- 
imation error as measured by n'' log^ n(/„/(/l"^' + /I^^' + fli^^ + 

fn ) — 1), with the /„ truncated after terms, for n = 50 . . 800 
(left) and for n = 776 . . 800 (right). 



a mere avatar of the product formula for the Gamma function, yields then easily 



3r(|) 



r(i + f)r(! + f)' 

The fluctuations of period 6 evidenced by Figure ^ are thus fully explained: one 



has 



J n ' J n ' J n 



o 



logn 



where the periodic function is (5R designates a real part) 



9r(f) 




n{n) =2(-l)" + 25R [uj- 



Expansions at z — rj, a primitive ith root of unity. Let 77 
Z = {1 — z/rj). The expansion of f{z) is now of the type 



exp(2i7r/£) and 



fiz) ^ /(yy) exp (Z + • • • + Z'-' \ogZ 



7I-1 



where Z^~^ logZ corresponds to the singular term in L\i{z^). Consequently, fluctu- 
ations start appearing at the level of terms of order in the asymptotic expansion 
of fn SiS n +00. The value of /(ry) is expressible in terms of Gamma values at 
algebraic points, as we have seen when determining f{uj). The coefficients in the ex- 
pansion also involve values of the V^-function {ip{s) = logr(s)) and its derivatives 
at rational points, which include C values as particular cases. 

Proposition 1. The probability that a permutation is made of cycles of distinct 
lengths admits a full asymptotic expansion of the form 



fn 



H 1 7T- (- log n - 1 - 7 + log 2) 

n 



(log' n + C3,i logn + 03,0) + 2(-l)" + 25R 



^ Prjn) 

r>4 



9r(i)a>- 



r(i+^)r(i-^) 



n' 
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with C3.1 and C3.0 as given by (jSSJ- There, Pr{n) is a polynomial of exact degree r — 
1 in \ogn with coefficients that are periodic functions of n with period D{r) — 
lcni(2, 3, . . . , r). 

Figure |21 displays the error of the approximation obtained by incorporating all 
terms till order included in the asymptotic expansion of /„. Fluctuations of 
period 12 (due to the additional presence oii = V— 1) start making an appearance. 



We now examine several combinatorial problems related to permutations, poly- 
nomials over finite fields, and trees that are amenable to the hybrid method. The 
detailed treatment of permutations with distinct cycle lengths can serve as a beacon 
for the analysis of similar infinite product generating functions, and accordingly our 
presentation of each example will be quite succinct. 

In the examples that follow, the function / whose coefficients are to be analysed is 
such that there is an increasing family of sets Z^^^ Z^'^\ . . . (ordered by inclusion, 
and with elements being roots of unity), attached to a collection of asymptotic 
expansions having smaller and smaller error terms. In that case, a full asymptotic 
expansion is available for the coefhcients of /. The general asymptotic shape of [z"]f 
involves standard terms of the form n^^ log'' n modulated by complex exponentials, 
since [z"]a{z/() — (~"'an- We formalize this notion by a definition. 

Definition 7. A sequence (fn) is said to admit a full asymptotic expansion with 
oscillating coefficients if it is of the form 



where the exponents increase to +cx) and each Pr(n) is a polynomial in logn 
whose coefficients are periodic functions of n. 

In a small way, Proposition^and the forthcoming statements, Propositionsl^HHl 
can be regarded as analogues, in the realm of functions of slow growth near the 
unit circle, of the Hardy-Ramanujan-Rademacher analysis (X, 2li ,20 of partition 
generating functions, the latter exhibiting a very fast growth (being of infinite 
order) as \z\ 1. 

4.1. Permutations admitting an m-th root. The problem of determining the 
number of permutations that are squares or equivalently "have a square root" is a 
classical one of combinatorial analysis: see Wilf's vivid account in ^40, §4.8]. The 
problem admits an obvious generalization. We shall let Umiz) be the exponential 
generating function of permutations that are mth powers or, if one prefers, admit 
an mth root. 

How many permutations have square roots? For the generating function 
112, we follow Wilf's account. Upon squaring a permutation r, each cycle of even 
length of T falls apart into two cycles of half the length, while an odd cycle gives rise 
to a cycle of the same length. Hence, if a has a square root, then the number of cy- 
cles it has of each even length must be even. By general principles of combinatorial 
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analysis, the exponential generating function of the number 112, n of permutations 
of n elements that have square roots satisfies 



7 

n2(z) y n2„^ = cosh(zV2)e" /3 cosh(zV4)e"^ 



n>0 



/5... 



(34) 



The series starts as 



exp 



1 , 1 + z 
— log 

2 ^ 1 - z 



cosh 



1 + 2 
1-2 



cosh 



m>l 

2m 



^2m 

2to 



n2(.) = i + z + |^ + 3^ + i2^ 



-.5 



its coefficients being EIS:A003483, with the quantity [z"]n2(2:) representing the 
probability that a permutation is a square. The GF of (|34|1 is given by Bender PI 
p. 510] who attributes it to Blum [S]. It is interesting to note that Bender mentions 
the following estimate from [S], 



(35) 



[^"]n2(^) 



and derives it by an application of the Tauberian theorem'* of Hardy, Littlewood, 
and Karamata. Accordingly, no error terms are available, given the nonconstructive 
character of classical Tauberian theory. We state: 

Proposition 2. The probability that a random permutation of size n has a square- 
root admits a full asymptotic expansion with oscillating coefficients. In particular, 
it satisfies 



(36) [2"]n2(2) 

where 



,G 



1 - 



logn Cg + (-1) 



4n 



2eG(-i)L»/2J ^ / logn 

^2 \^ j^5/2 



(37) 



C2 



n-»(i) 

k>l ^ ' 

V ( -i- -tanh(-i- 



k>\ 



= 1.22177 9515192536 

C3 = -12 + 161og2 + 47 + 2c2. 



Proof. First a rough analysis suffices to see that 



and 



^2iz) 



1 + 2 



z^l V 1 - 2^ 2^-1 V 2 

In order to refine these expansions, introduce the normalized tangent numbers by 
tan 2 = rm2^'"^*, so that logcos(2) = Tm~i . 



^Regarding Tauberian side conditions, Bona, McLennan, and White ^ prove by elementary 
combinatorial arguments that the sequence 112, n is monotonically nonincreasing in n. 
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The usual exp-log reorganization of the 112 series yields 



' \m>l 

In passing, this provides for G, the fast convergent series 

m>l 

on which the numerical estimate of (|37|l is based. 
Next, take out the e*^ factor, leading to 



(38) n^iz) = e«yi±| exp ^"^^ y^T^™-! [Li^™ (z'^) - C(2m)] 

At z = 1, the largest singular term in the exponential arises from Li2(z'*), the 
contributions from the other polylogarithms being of smaller order 

U^iz) ^ e^'J-^ (1 + 2(1 - z) log(l - z) + c*(l - z) + O ((1 - z)^ log2(l - z))) 

where, here and later, c* designates a computable constant that we leave unspecified 
for the purpose of readability. Similarly, at z = — 1, we find 



n,{z) e«^ii^ (1 + O ((1 + z) log(l + z))) . 
At z = i = \/— T (hence at z = — i, by conjugacy), we have 

Tl2{z) ^ (P*(l - z/t) + 2(1 - z/i) log(l - zA) + 0((1 - z/z)3/2)) , 

with P* an unspecified polynomial that does not leave a trace in the coefficients' 
expansion. 

The singular contribution to 112 arising from any root of unity C of order > 8 is 
at most (0(1 — z/()^ log(l ~ z/()), which translates to an 0(n^^) term. The proof 
of H3t)|) is then completed upon making c* in the expansion at 1 explicit, which 
introduces the new constants 02,03. Existence of the full expansions finally follows 
from the usual analysis of polylogarithms of powers, taken at roots of unity. □ 

How many permutations have an mth root? Like the previous one, this 
problem is briefiy mentioned in Bender's survey [3]. We follow again Wilf's expo- 
sition |4(J[ §4.8]. For a pair £,m of positive integers, we define ((£, to)) to be 

((^, to)) := lim gcd(F,TO). 

(Thus, ((£, to)) gathers from the prime decomposition of to all the factors that 
involve a prime divisor of £.) The characterization of permutations that are mth 
powers then generalizes |S| EH El EDI • a permutation has an mth root if and only 
if, for each ^ = 1, 2, . . ., it is true that the number of cycles of length £ is a multiple 
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of {{£, m)). This observation leads to an expression for the corresponding generating 
function. Indeed, define the "sectioned exponential" , 

n>0 ^ ' j=0 

SO that expi{z) — and exp2(2;) = cosh(z). The exponential generating function 
of permutations that are mth powers is then 

(39) n„(z) = H'^^Pa^,™)) T 



The generating function of H39() has been investigated by Pouyanne whose 
paper provides the first order asymptotic estimate of [z"]n„i. There is a funda- 
mental factorization, 

(40) n™(z) = A„.(z)-B™(z), 

where A„x gathers from the product (|39|l the numbers t that are relatively prime 
to m and B{z) gathers the rest. 

The factor A^, is found by series rearrangements to be an algebraic function 
expressible by radicals, 

A^{z)= [] (l_z'=)-^w/^ 

k I m 

with /i(fc) the Mobius function. For m = 2 (square permutations), this is the 
ubiquitous prefactor ^2(2) = y^(l + z)/(l — z). For m = 6, for instance, the 
prefactor becomes 



The factor B^a is a transcendental function that admits the unit circle as a natu- 
ral boundary |32| . It is expressible as an infinite product of sectioned exponentials: 



(41) -B,„(z) = W exp((^ ( y 



Ecd(*,m)>l 



For m = 2, this is the infinite product of hyperbolic cosines. For m = 6, one has 

B^{z)^ n ^^Pefy) n exp^fy) n ^^psfy 



l=Q mod 6 ^ ' 1=2,4 mod 6 ^ ^ £=3 mod 6 ^ 

These singular factors can be analysed just like in the case of 112 by an exp-log 
transformation. One first observes that the limit value 5^(1) is well defined, since 
the infinite product converges at least as fast as ^(1 + ^~^)- It is seen next that 
singularities are at roots of unity, and the radial expansions can be computed in 
the usual way from the polylogarithmic expansion. We can now state a (somewhat 
minor) improvement over |32j : 
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Proposition 3. The probability that a random permutation of size n has an mth 
root admits a full asymptotic expansion with oscillating coefficients. To first as- 
ymptotic order, it satisfies 

L J ™v ; ^i-<p(m)/m' " r{ip{m,)/m^ 

where ip{m) is the Euler totient function and 
(42) i3™(l)= Tl expp,„)) - 



k I m 



cd(*,m)>l 



In particular, when m — p is a prime number, one has 

(43) [^"]n,(.) ^ := ^.^J^ n -P. [t. 



For instance, we find 

nn6(z) = S6(i 



12 ^_2/3 ^ 2J^^ (^p-.f (n+D) „-5/6 + 0^1 



.r(i) r(i) 

and 143|l improves upon an early estimate of Turan ^| Tfi. IV] . 

4.2. Pairs of permutations having the same cycle type. Given a permutation 
CT, its cycle type is the (unordered) multiset formed with the lengths of the cycles 
entering its decomposition into cycles. For a permutation of size n, this type can 
be equivalently represented by a partition of the integer n: for instance (2^,5,7^) 
represents the profile of any permutation of size 25 that has three cycles of length 2, 
one cycle of length 5, and two cycles of length 7. The probability that a random 
permutation of size n has profile 

^ ^ (•-^ni^2"^ . . .), with lrii + 2n2H = n 

is, by virtue of a well-known formula 8, p. 233], 

Z>1 

corresponding to the generating function in infinitely many variables 



(44) <^{z\ xi,X2, • ■ •) = exp yxiz + X2 — + x^, — + ■ ■ 

which is such that — [z'^x^^x^'^ ■ ■ •]<!>. 

In this subsection, we estimate the probability that two permutations of size n 
taken uniformly and independently at random have the same cycle type. (Each pair 
of permutations is taken with probability 1/n!^.) The quantity to be estimated is 
thus 

Wn = ^ (a^) , 

7r h n 

where the summation ranges over all partitions tt — (1"S 2"^, . . .) of n. This prob- 
lem was suggested to the authors by the reading of a short note of Wilf , who 
estimated the probability that two permutations have the same number of cycles 
(the answer to the latter question turns out to be asymptotic to l/(2-^7r logn)). 
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Given ijl^ . it is not hard to find the generating function of the sequence (Wn)- 

°° / k\ n 

(45) m-):=E^"^" = n^ p -here /(z) = ^ _, 

n>0 fe=l ^ ^ n>0 

the reason being that 

^w = EE n72^ = ni E 



>1 \ni>0 



The function written /(z), which is obviously entire, is a variant of the Bessel 
function Iq (see, e.g., |31]): I{z) — Io{2y^). Also, from H45() . the expansion of 
H^(z) is readily computed: one has 

m-) = l + - + 2|^ + 14^ + 146^ + 2602|^ + ..., 

where the coefficients are EIS:A087132 ("sum of the squares of the sizes of con- 
jugacy classes in the symmetric group S„"). 

Proposition 4. The probability Wn that two permutations of size n have the same 
cycle type satisfies 



k>l 



= 4.26340 35141 52669. 



Furthermore, this probability admits a full asymptotic expansion with oscillating 
coefficients. 



Proof. We shall only sketch the analysis of the dominant asymptotic term, the rest 
being by now routine. From the expression of W{z), the exp-log transformation 
yields 



> fc>i 



W{z) = exp I E ( ^ ) ) where H{z) := log (/(z)) . 



Let he = [z^]H{z), that is, H{z) = hiz^ = z — jz^ + ■ ■ ■ (the sequence {h(i\^) is 
EIS:A002190, which occurs in the enumeration of certain pairs of permutations 
by Carlitz). Then 



W{z) = exp E E pi = E heU2e{z'^) 

\fc>l^>l J \£>1 

This expression ensures that the hybrid method can be applied at any order, with 
W{z) being of global order a = 0. The first term of the asymptotic estimate is 
provided by the factorization relative to Z = {1,-1}, namely W{z) — P{z)Q{z), 
with 

P{z) = exp (Li2{z) ~ jLiiiz"^)! , Q{z) = exp | hiLhij 
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Since Q{z) is clearly on the closed unit disk, the hybrid method applies: with 
the notations of Theorem [21 one can take s = 4, to the effect that uo — 2. Using 
the algorithmic scheme of Section 13 we find 

Wiz) = W{1) [l - (1 - z)L + 1(1 - z^L + i(l - z^L' + 0((1 - zf) 

W{z) = W{-l)[l + 0{{l + zf)], 

2:^ — 1 + 

with L = L{z) = log(l — z)^^. Theorem 13 then directly yields W„ ~ and 
further asymptotic terms can easily be extracted. □ 

Regarding other statistics on pairs of permutations, it is well worth mentioning 
Dixon's recent study of the probability that two randomly chosen permutations 
generate a transitive group For the symmetric group S„, this probability is 
found to be asymptotic to 

^11 4 23 

n n? rfi 

and, up to exponentially smaller order terms, this expansion also gives the proba- 
bility that two random permutations generate the whole symmetric group. In that 
case, the analytic engine is Bender's theory of coefficient extraction in divergent 
series 0. 

4.3. Factorizations of polynomials over finite fields. Factoring polynomials 
in Q\X] is a problem of interest in symbolic computation, as it has implications 
in the determination of partial fraction expansions and symbolic integration, for 
instance. Most of the existing algorithms proceed by a reduction to a few factoriza- 
tions of polynomials with coefficients in some finite field (with q a prime power): 
see the books by Berlekamp, Knuth, or von zur Gathen and Gerhard [H 1381 1^ . 
Several of the algorithms for factoring polynomials over finite fields involve what 
is known as the distinct degree factorization. To wit, let (j) £ '^q[X] be completely 
factored as 



' 
] 

where the Uj are distinct irreducible polynomials. The distinct degree factorization 
produces the decomposition of the square-free part of (j) under the form 

squarefree(0) = f3if32 ■ ■ ■ (3n, where /3fe — tj . 

deg(t3) = fc 

We refer to the paper jj^j by three of us for details. 

The distinct degree factorization of as defined above coincides with the com- 
plete factorization of (j) if and only (j) has all its irreducible factors of different 
degrees. For the purpose of analysis of algorithms, it is then of interest to quantify 
the probability of that event. This problem has been considered independently by 
the authors of jl2[ll3) and by Knopfmacher & Warlimont in 1241 . 

First, without loss of generality, we restrict attention to monic polynomials. If q 
is the cardinality of the base field Fg, the number of monic polynomials of degree 
n is P„ = g", with ordinary generating function 

P(z) = ^g"z" = -^. 
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The unique factorization property impUes that the class V of aU (nionic) polyno- 
mials is isomorphic to the class obtained by applying the Multiset construction 
to the class X of (nionic) irreducible polynomials. This, by well known principles 
of combinatorial analysis 0] I15| translates to an equation binding the generating 
function I{z) of irreducible polynomials, namely, 

(47) P{z) = n (1 - z^-'" = cxp (l{z) + \l{z^) + \l{z') + •••), 

n— 1 ^ ' 

where /„ is the number of monic irreducible polynomials of degree n. The solution 
is obtained by taking a logarithm of the second form and applying Mobius inversion: 



(48) /(,)^^^log 



k 1 — qz^ ' 
k>i ^ 

and, as a consequence, the number /„ admits the explicit formula (already known 
to GauB, see references in ^1 p. 46]): 

(49) In^-Y] mi"^^ = — +Rn where |i?„| < q^'\ 

n ^ — ' n 

k I n 

In passing this approach also shows that 

/(z) = log^^ + i?(z), 

where R{z) is analytic in \z\ < q^^^^, so that I{z) has a logarithmic singularity at 
z = q~^. 

Let Dn be the number of (monic) polynomials composed of irreducible factors 
of distinct degrees. The corresponding generating function, D{z) satisfies 

oo 

(50) D{z) = n (1 + ^"^") ' 

n=l 

and the probability of the event over polynomials of degree n is 

It is in fact the simplification of this problem for fields of large cardinalities that 
originally led Greene and Knuth to consider the problem of enumerating permuta- 
tions with distinct cycle lengths (see (T31 for context). We have: 

Proposition 5. The probability that a random polynomial of degree n over ¥q has 
all its irreducible factors of distinct degrees satisfies the asymptotic estimate 

(51) £^^5iq) + o(- 
where the constant S{q) is given by 

fe>l ^ ^ ^ 



k\Ik 



In addition, this probability admits a full asymptotic expansion with oscillating co- 
efficients. 
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Proof. We start with a rough analysis based on comparing D{z) and P{z). The 
second form of H47|) imphes 

(53) ^=n(i+^»'?""-")(i 

^ n>l 

At z = 1, the general factor in the product satisfies, for large n, 
(54) 

(l + /„q-")(l-r")'" = (^l + i + 0(g-"/^))exp(^-i + 0(g-"/2 

- l + i + 0(n-3), 

given the known form 149|) of /„. This is enough to ensure the convergence of the 
infinite product defining d{q) and proves that 

which is evidently compatible with the asserted estimate (|51() . Also, the esti- 
mate H54|l points to the fact that convergence of the infinite product giving i5(g) is 
very slow, but also suggests the modified scheme (compare with H24|) '). 

(55) 6{q) = e--' n (^Tq^(i - ^-'y-'^'^' 

which now exhibits geometric convergence. 

In order to complete the derivation of H51|l , we apply the exp-log transformation 
to ISOI): 

(56) D (-] = exp I ^ IziZl^A„(z) j , A„(z) ^ (/„(z""z")'" . 

\m>l J n>l 

Also, we have Ai{z) — I{z/q). Taking out the factor e'^^'^^-' — e^'^^^'^'^ gives 
15(f) =e^(^/').exp| ^H^A„(z) 

-T^expfEt^log^3^]expfE 

— — — ^"(^) 

I m>2 



^ ^ ■A„(z) 



TO 
2 



1 — Z \ ^ — ' TO 



where A{z) is analytic in \z\ < q^^^. 

Next, we observe from (I49II that 

A„,(z) =Li„,(z™) + S'™(z), 

where Smiz) is analytic in \z\ < q^l"^. The existence of a factorization fulfilling 
the conditions of Theorem 13 now follows from the usual split of the sum in the 
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exponential of (|57|l . last line: 

> ^ A,„ z = > '- A„. z + > '- A„ z . 

m>2 m=2 m,>M 

The first split sum is singular at roots of unity of order < M and admits a log-power 
expansion of type O* for any t; the second split sum has its nth coefficient that 
is of order 0{n~^~^), and is accordingly C^~^. By Theorem |21 these consider- 
ations imply that [z"]D{z/q) admits a full asymptotic expansion with oscillating 
coefficients. 

It is not hard then to track the order of the term that corrects the dominant 
asymptotic regime: in essence, it arises from the factor 

1 C-il^h(z') 



1 

whose coefficients are 0(ri~^). □ 

A variant of this problem is also considered in |13[ : determine the probabil- 
ity that a random polynomial of degree n is such that its square-free factorization 
consists of polynomials of distinct degrees. The generating function is then 

n>l ^ 

By devices entirely similar to those employed in the proof of Proposition |31 one 
finds that the corresponding probability satisfies 

%^%) + o(i 

where the constant 6(q) is 

^('^) - n (i + ^) (i-^-^r^ - n {^^TTw^^' - '^''y-^'^') ■ 

The treatment of is based on a crude application of the hybrid method, that 
of |24| relies on elementary bounds and coefficient manipulations, but it is not 
clearly applicable to derive asymptotic expansions beyond the first term. 

4.4. Dissimilar forests. Fix a class T of rooted trees. An unordered forest is a 
multiset of trees in T, and we let T be the class of all forests. A forest is said to 
be dissimilar if all the trees that it contains are of different sizes. We let £ denote 
the class of dissimilar forests. 

For instance, the collection T of all plane trees has generating function 



(58) T{z) = ; 

a well-known fact jl5[ I17L I25L I35L QUI ■ The generating function of unordered forests 
is given by a construction analogous to l|Tf|) : 



F{z) 



= n(i - -pfE^^(-^-)) 

n>l \k=l I 

= l + z + 2z^ + 4:Z^ + lOz^ + 26z^ + 77 z^ + 235z'^ 



HYBRID METHOD 



27 



(the coefficients constitute EIS:A052854). The generating function of dissimilar 
forests, which is our actual object of study, is 

oo 

E{z) = n(l + T„z") 

n=l 

= l + z + z^ + 3z^ + Tz-* + 21z^ + 63z6 + 2082'^ + • • • 

(currently unlisted in |34|). 

The coefficients [z"]T(z) are the famous Catalan numbers of combinatorial the- 
ory, 

1 /2n - 2 



T 

n \ n — I 

the asymptotic estimate being in agreement with the square-root singularity that 
is visible in T(z). An easy analysis of F{z) then shows that, as z ^ j, 

1 



F{z) = K - -K\/l -Az + --- 

(59) 



K := exp I > — (1- VI -41-'=) )= 1.7160305349 22281, 
resulting in the estimate 

4"-i 

F, - K^= 



Proposition 6. The probability that a random forest of rooted plane trees having 
total size n is dissimilar is asymptotic to the constant 

(60) K = — with L:= y ^ ^ > - 



K ^-^ m ^-^ \ n\ ri — 1 

m=2 n>l ^ ^ 

This probability admits a complete asymptotic expansion with oscillating coeffi- 
cients. 



Proof. We limit the calculation to first order asymptotics. As usual, start from the 
exp-log trick, which yields 




(61) Eiz) = exp > ^ ^ A„,(z) 



m=2 



where the functions Am that play the role of special polylogarithms are defined by 



(62) A„(z):=5:(i^2" 2 



n>l 



n\n — l 



z 



Observe that the functions Am(z) are A-analytic, by virtue of a theorem asserting 
the closure of functions of singularity analysis class under Hadamard products |10l 
ITS] . (The function A,„(z) is roughly comparable to Li3„/2(4'"z'")-) Splitting the 
sum in (|61|l according to m < M and m > M, for an arbitrarily large M, yields a 
factorization of E{z) that fulfills the conditions of Theorem[21 Hence the existence 
of a full expansion for coefficients is granted. 
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Generating function 



n>l 



n>l ^ 
n>l ^ 



g-oo = n 1 + 



Type 



Partitions dist. summ. 



Perm. dist. cycle lengths 
Forest dist. comp. sizes 
Perm, same cycle type 



Set part. dist. block lengths 



Coefficients 



oc n ^''^62 



oc 1 



-3/2 



Mixed regimes |23| . 



Figure 3. Some typical infinite-product generating functions, 
their "types" , and the asymptotic shape of the corresponding co- 
efficients. 

For first order asymptotics, it suffices to analyse what happens as z — > (^)^. We 
find, from H()f|) and with the notations of (|60|l 

f 



E{z) = ( f 



which implies the stated estimate. 



-VI -4z + 0(f -4z) 



□ 



The analysis obviously extends to any simple family of trees, in the sense of Meir 
and Moon |15[l29j . including the case of Cay ley trees that are enumerated by 

5. Conclusion 

As demonstrated by our foregoing examples, several infinite-product generating 
functions occurring in combinatorics are amenable to the hybrid method, despite 
the fact that they admit the unit circle as a natural boundary. For comparison 
purposes. Figure 13 displays a few typical infinite-product generating functions and 
the general asymptotic shapes^ of their coefficients. 

Perhaps the most well-known generating function in this range is 



whose coefficients enumerate partitions into distinct summands. In such a case, the 
analysis of coefficients, first performed by Hardy and Ramanujan, is carried out by 
mean of the saddle point method, in accordance with the fact that the function 
Qq{z) gets exponentially large at its "main" singularity z = 1. This and similar 
cases are best treated by means of Meinardus' method jlj. 
The function 



ri>l ^ ^ 



We use 'oc' to represents asymptotic proportionality, that is, asymptotic equivalence up to an 
unspecified multiplicative constant. 
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is also fast growing near z = 1 and is best transformed by means of convergence 
factors: 

This expression shows a growth similar to that of the function exp((l — to 
which the saddle point method is applicable. All computations done, we find 

g_i/2,„'-G^expg^) with G=n + 
The next three examples of Figure 01 are relative to functions of the form 

n>l 

with a < — 1. These are typical cases where the hybrid method is applicable. The 
function Q-i is exactly the one encountered when enumerating permutations with 
distinct cycle lengths. The function Q-3/2 provides a simplified analytic model of 
dissimilar forests, while Q-2 behaves like the generating function of permutations 
having the same cycle type. The general pattern is, for a < — 1: 

At the other end of the spectrum, we find the function 

n>l ^ ' ^ 

which is the exponential generating function of set partitions, all of whose blocks 
are of distinct sizes. Clearly, this is an entire function to which the hybrid method 
is not applicable. Coefficients have been precisely studied by the seven authors 
of with some (but not all) of the regimes being accessible to the saddle point 
method. 

Finally, Gourdon has developed in [18j methods leading to asymptotic expansions 
of a shape similar to the ones of the present paper in relation to a refinement 
of Golomb's problem posed by Knuth in [551 Ex. 1.3.3.23]. He showed that the 
expected length of the longest cycle in a permutation of size n admits an asymptotic 
representation which starts as (a; = e^"/^) 

A" + - ^ + (|I - U-^r)w + (i£ + i(-ir + i-^-" + X"^);^, 

with higher order roots of unity appearing at higher asymptotic orders. The proof 
has analogies to the techniques of the present paper, with the additional need of 
precise expansions that describe the behaviour of the truncated logarithmic series 
on the unit disc. 

Acknowledgements. Bruno Salvy's equivalent program has proved invaluable 
in permitting us to make several of our asymptotic expansions explicit. 
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